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Abstract 

An improved efficient method of calculating the hard bremsstrahlung correction to e~^e~ — > 4/ 
for non-zero fermion masses is presented. The non-vanishing fermion masses allow us to perform 
the phase space integrations to the very collinear limit. We therefore can calculate cross sections 
independent of angular cuts. Such calculations are important for background studies. Results are 
presented for the total and some differential cross sections for e~^e~ — > ud^~i?^ and the corresponding 
hard bremsstrahlung process. The latter is of particular interest for a detailed investigation of the 
effects of final state radiation. In principle, the process e~'"e~ ud^~i>^^ is also interesting since 
it helps to set bounds on possible anomalous triple and quartic gauge boson couplings involving 
photons. The size of mass effects is illustrated by comparing the final states ud/i" 1/^(7), cs^~9^j_{^) 
and udT~UT-{'^). 
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1 Introduction 



Precision measurements of properties of the intermediate gauge bosons Z and W have deep- 
ened our understanding of the electroweak interactions and consohdated the vahdity of the 
electroweak Standard Model (SM) considerably in the past decade. After the convincing suc- 
cess of LEPl and SLD experiments in pinning down the properties of the Z-resonance we 
expect further advances in measurements of the properties of the W boson which are not yet 
known with comparable precision. As the SM very precisely predicts the mass and the width 
of the W a high accuracy determination of these parameters is of crucial importance, because 
they allow us to improve indirect bounds on the Higgs mass or on new physics beyond the SM. 
The precision of ongoing measurements, single-VT production at the hadron collider TEVA- 
TRON, and VT-pair production at LEP2 are limited by statistics and/or by lack of detailed 
theoretical understanding. 

The proper analysis of VT^-pair production at LEP2 and later at future high energy e^e~ 
linear colliders requires the accurate knowledge of the SM predictions including all rele- 
vant radiative corrections. What we need is a detailed understanding of the production 
and subsequent decay of W pairs, including the background processes and photon radia- 
tion: e^e~ 4/, 4/7, 4/77, ■ ■ -, where 4/ denotes a possible four fermion final state. The 
lowest order theoretical results for all the possible four-fermion final states have been already 
implemented in several Monte Carlo event generators and semi-analytic programs, which have 
been thoroughly compared in [Q]. Most of the programs include some classes of radiative cor- 
rections such as the initial- and final-state radiation. Coulomb corrections, running of the fine 
structure constant, etc. While presently available e^e" 4/, 4/7 matrix elements are precise 
enough for the analysis of LEP2 data [Q], at future linear colliders, a much better knowledge 
of the radiative corrections will be necessary because of the high statistics expected at these 
accelerators and because radiative corrections get more significant at higher energies. 

The complete one-loop electroweak radiative corrections to the on-shell pair production 
including soft bremsstrahlung were calculated in 0. The hard bremsstrahlung process e+e~ 
W~^W~'j was included in and ||^. For the process e~^e~ W^W~ 4/ of actual interest 
to the experiments only partial results are available. We refer to for a recent review of the 
status of precision calculations for this case. 

Sufficiently above the pair production threshold, for most of the present applications, it 
seems to be sufficient to take into account corrections to the double-resonant diagrams only, 
i.e., e~^e~ 4/ via virtual W^W~ intermediate states. The validity of this approximation 
has to be controlled by more complete calculations, however. From a theoretical point of view 
it is certainly necessary to evaluate the complete 0{a) radiative corrections for the different 
channels of the 2 — >■ 4 fermion reactions. However, despite of the fact that some progress in 
calculating the complete virtual one-loop electroweak radiative corrections to e^e~ —>■ ud^'v^ 
have been reported in 0, the final result of such a calculation is still missing. Concerning 
the real photon emission, the situation looks much better. The hard bremsstrahlung for 
four-fermion reactions mediated by two resonant W bosons was calculated in ^j. A similar 
calculation, extended by an inclusion of collinear effects, was presented as a package WWF 0. 
The complete lowest order result for e+e~ — >• e~Ueud'y was presented in [jlO| and calculations 



of e^e Af'-f for an arbitrary final state were reported in [^. Results on bremsstrahlung 
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for purely leptonic reactions have been published in and most recently predictions for all 
processes e~^e~ — 4/7 with massless fermions have been presented in [|12| . 



At a future linear collider, the proper treatment of the collinear photons will be crucial and 
it requires to take into account the fermion masses appropriately. Therefore, in the present 
paper, we propose an efficient method of calculating the hard photon bremsstrahlung for 
four fermion production in e^e~ annihilation without neglecting the fermion masses. The 
phase space integration can therefore be performed to the very collinear limit. This allows for 
calculating cross sections independent of angular cuts and estimating background contributions 
coming from undetected hard photons. We present results for the total and a few differential 
cross sections for the channel e~^e~ —>■ udfi~Un and the corresponding bremsstrahlung process. 
The latter is particularly suited for a detailed investigation of effects related to final state 
photon emission, since the muons appear well separated from photons in the detectors. In 
particular it seems to be interesting to study the infiuence of final state radiation on the W mass 
measurement via this channel. Having the final state photon resolution in e~^e~ — > udfi~Ufj_'j 
could also make it possible to investigate the quartic jVWW couplings (V = 7,2'), which 
are absent on the Born level of 4/ production. Of course, besides the new quartic couplings 
there are additional triple ''yWW vertices as well. In the soft photon limit, we can perform 
the integration over the soft photon phase space analytically and demonstrate the cut-off 
independence of the combined soft and hard photon bremsstrahlung cross section. We finally 
will illustrate the importance of mass effects by comparing the channels where ud is replaced 
by cs. Similarly, we may replace the /i by a r lepton. 



2 Method of calculation 



In this section, we present a method for calculating the matrix elements of a two-fermion to 
four-fermion reaction and an associated bremsstrahlung photon. The method is an extension 
of the helicity amplitude method introduced in [^] to final states of arbitrary spin. 

As in 1^, we use the Weyl representation for fermions where the Dirac matrices 7'^,/i = 
0,1,2,3, are given in terms of the unit 2x2 matrix I and Pauli matrices ai,i = 1,2,3, by 

with a!^ = (/, ±(Tj). In representation (|I]), the matrix 75 = 'i7'^7-'^7^7^ and the chiral projectors 
P± = (1 ±75)/2 read 

-^(V?)^ ^-^(il)^ --(S?)' (^) 

A contraction of any four- vector a'^ with the 7^^ matrices of ([I|) has the form 

The 2x2 matrices can be expressed in terms of the components of the four-vector a'^ by 
+ _ ( a° — —a^ + ia'^ \ _ _ ( + — ia? \ . . 



In representation (|1]), the helicity spinors for a particle, ^(p, A), and an antiparticle, v{p, A), 
of four- momentum {E, p) and helicity A/2 = ±1/2 are given by 





A 


Pl 


'e- 


A 


Pl 



and the helicity eigenstates ^(p, A) can be expressed in terms of the spherical angles 6 and 
of the momentum vector p as0 

^^P'+^) = [ e^'^ sin 9/2 ) ' ^^P' = [ cose/2 ) ' 

For simplicity we use real polarization vectors which are defined again in terms of 9 and 

e^(p, 1) = (0, cos 9 cos 0, cos 9 sin 0, — sin 6^) , £^^(p, 2) = (0, — sin 0, cos 0, 0) (7) 
e^(p, 3) = 7 (/?, sin 9 cos 0, sin 9 sin 0, cos 6^) , (8) 

where the longitudinal polarization component of (|]) is defined exclusively for a massive vector 
particle of energy m7 and momentum m7/3. We could use complex polarization vectors in the 
helicity basis as well, if we were interested in definite helicity polarizations. 

A polarized matrix element is calculated for a given set of external particle momenta in a 
fixed reference frame, e.g. the center of mass system (c.m.s.) of the initial particles, where 
the initial momenta are parallel to the z axis. 

Fermion masses are kept nonzero. The mass effects play an essential role in the bremsstrahlung 
reactions whenever a collinear photon is emitted. They also are important for tree level 
reactions with identical particles in the initial and final state, where a photon exchanged in 
the t-channel approaches its mass-shell. Moreover, by keeping the fermion masses finite the 
Higgs boson exchange can be incorporated in a consistent way. 

In order to speed up numerical computation, we decompose the Feynman graphs into factors 
which depend on a single uncontracted vector-index and therefore may be considered as gen- 
eralized polarization vectors. They can be easily computed and used as building blocks of 
other graphs. For example, the coupling of an internal gauge boson to the external fermions 
may be considered as the generalized polarization vector Sy which is defined as 

e^y ip,,p2, Ai, A2) = Ai)7. {g^^^P- + 9^v ^P+) MP2, A2) 

+ (mi^{,+) - mjy^) )] / - Ml) (9) 

where '0i(pi, Ai) = (-^l, ^[^) and ip2{P2, ^2) = ("025 i^i^) spinors for a particle or an antipar- 
ticle of four-momentum pi, mass rrij and helicity Aj, as defined in (^. We have denoted the 
chiral couplings of the '0i'02^ vertex by gy^\ Dy'{q) is the photon propagator in the Feynman 
gauge or the massive gauge boson propagator in the unitary gauge and q = ±pi ± p2 is the 

^Note that our phase convention differs from the one chosen in [Q. 
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four-momentum transfer. The +(— )-sign corresponds to an outgoing (incoming) particle. In 
the case of a photon propagator we have My = and only the first term in the square brackets 
on the right hand side of (^ is present. 

The photon emission from any of the external fermion legs of the ijJiijj2V vertex can be taken 
into account by defining two other generalized polarization vectors: 

£'^v (Pi,P2, k, Ai, A2, A) = Df^^iq + k) ipiipi, Xi)g^i^{k, A) ^ f\t 2 

[±Pi + ftj — nil 

= ^^^r(9 + k) [g^^^i^i (2pi ■ e t k^e') a+^f + g^y^i' {2p, ■ e t k~e^) a'^i] (10) 

where the upper (lower) sign is assumed if ipi represents an outgoing particle (incoming an- 
tiparticle) and 

5^, (pi,P2, k, Ai, A2, A) = D^\q + k) Mpi, Ai)7. {g^F^P- + gt^^P+) 



D!i.^{q + k) WHiat (~2p2 ■ e ± A;-£+) + g^yHi' a' (-2p2 ■ e ± k+e-) ^ 



2p2 ■ k 
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where the upper (lower) sign has to be taken when ip2 represents an incoming particle (outgoing 
antiparticle). In (|10]) and ([TT|) , e^{k, A) is the photon polarization vector as defined in (|^) and 
g^i are the photon couplings to ^pi. 

If we contract the triple gauge boson coupling 

^^WWVM^P2,P3) = 9WWV [{Pl - P2) + {P2 - PsTg"' + (P3 " J>l) VI , (12) 

where pi,P2 and p^ are the incoming momenta of the W^l, W~ and the neutral gauge boson 
V^, = 7, Z", respectively, with two (generalized) polarization vectors, say 82, and with a 
gauge boson propagator we will obtain another generalized polarization vector, e.g. 

£^(l,2) = D^^(g)r('J,^^)ere^. (13) 

In (pr3|), 1 and 2 stay for the four-momenta and polarizations. With the help of the generalized 
polarization vectors (P-p|) and (|13D, the amplitude corresponding to any Feynman diagram 
of a process e+e~ 4/(7) may be represented by one of the scalar functions F2n+i, -E3 or 
we are going to define now. A fermion line containing n + 1 couplings to gauge bosons and 
n fermion propagators sandwiched between external spinors can be represented by the scalar 
function 

F2n+i (^ipi,g['^^ei,g['^ei,pi,Pi,mi,g^2^^6j,gi'^62,pt,P2,m2, 

{A-'p- + Pi-'p.) {A-'p- + p.) ■ ■ ■ 

^ ^ PI — mf ^ ^ P2 — m2 

^4+1 {gitP- + 9i^>^) ^2, (14) 

5 



X 2 

Pi - m: 



where we have suppressed polarization indices. In the representation the algebra of 4 x 4 

matrices in (|l^ can be easily reduced to the algebra of 2 x 2 matrices. This kind of reduction 
has been used already in (^|-|n|). Utilizing the 2x2 algebraic representation speeds up the 
numerical computation and allows for a simultaneous calculation of the 7 and Z contributions. 
We have indicated the use of the reduced form by writing explicitly the dependence on the 
2x2 matrices ef,pf, defined according to (^) , on the left hand side of ([T^) . The general form 
of this reduction is quite a complicated formula. Therefore we restrict ourselves to present an 
example of the function F3 which can be written as 

+9[''^9i''^i'ietplet^i' + gtWi^x e\pXe~^i,{ + m^g^^^^^i' e^et^l' . (15) 

A contraction of the triple gauge boson coupling (|l^) with three polarization vectors £1,^2 
and £0 can be considered scalar function 

^3 [Pl,£l,P2,£^2,P3,£^3) = 

9wwv{{p\ -p2)-£3 £i-£2 + {P2 -P3)-£i £2-£3 + iP3-Pi)-£2 Si-^s] ■ (16) 
Similarly, the quartic gauge boson coupling 

rflf = gv.v.ww {g^'g"" + g^'^g"" - ig^'gn , (17) 

where the vector indices /z, z/ are associated with the neutral gauge bosons Vi and V2 and 
p, cr with and W~ , if contracted with four polarization vectors 61,62,6^ and £4 can be 
considered as another scalar function 

Ei {£i,£2, £3, £i) = gviV2WW (^i"^3 £2- £4 + £\-£^ £2- £3 — 2ei ■£2 £^3 "^4) • (18) 

Note that similar scalar functions can be defined for the Higgs boson by replacing the vector 
boson coupling and propagator in (|9HTll) and ([T3|) by the Higgs couphng and propagator. 

Functions (|T^), (|16D and (0) can easily be implemented in a FORTRAN program and computed 
numerically for any specific set of the particle momenta and polarizations. The Fortran 90 
language standard which contains a number of new features and intrinsic functions especially 
for array manipulations is particularly suitable for this task. 

The method described above can be used to calculate the matrix element of any process of 
e^e~ annihilation into four fermions and a photon. Actually, the method is quite general 
and can be applied to any 1 n tree level reaction, not necessarily in the framework of the 
standard model, with massive or massless fermions and/or bosons in the final state. However, 
it may happen that one or a few extra functions will have to be defined in addition to those 
defined in (pi^D, (|TCp and (|TBp. Practically the only limitation of the method is the feasibility of 
phase space integration which is performed numerically by applying the Monte Carlo method. 

3 Application to e^e~ ud/i^u^ and e+e~ udfi~Un^ 

Let us demonstrate how the method of Sect. ^ works in case of the tree level four-fermion 
reaction 

e^iPi, Ai) + e"(p2, A2) u{p3, A3) + J(p4, A4) + /i~(p5, A5) + i^^^iP6, Ae), (19) 
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where the four-momenta and hehcities are indicated in parenthesis. The Feynman diagrams 
of the process (|l^) are shown in Fig. [l|. Although the fermion mass effects are irrelevant for 
the reaction (|T^), we keep masses finite for the sake of illustration. However, we neglect the 
Higgs boson exchange contribution which is suppressed by ratios of the fermion masses to the 
W boson mass. 




(6),(7) 




d (8), (9) 

Figure 1: The Feynman diagrams of reaction 




We define the necessary generalized polarization vectors using (H) 



£^7(Pi,P2, Ai, Aa) 
£^z(Pi>P2, Ai, Aa) 

£^tV+(P3,P4, As, A4) 
^(^~(P5,P6, A5, Aq) 



^^''(^12) VliPl, >^l)lug'yeU2{P2, A2), 

^r(Pl2) MPU Ai)7. {g^'JP- + g'^ZeP^) MP2, A2), 

DwiPsi) usips, hhugwP-ViiPi, A4), 
DwiP5e) U5{p5, X^)iugwP-V&{.PQ^ Ag), 



(20) 



where we have introduced the shorthand notation = pi +P2, P34 = Ps +Pi and psg = p^ +Pq', 
g-ye^g'^Ze ^"^^ gw are the standard model couplings. We use constant widths for the massive 
gauge bosons. They are introduced through the complex mass parameters My = rriy — imv^ 



V 



in the propagators Dy ,V = W, Z. However, we keep a real value of the electroweak mixing 
parameter sin^vy- This simple prescription preserves the electromagnetic gauge invariance, 
also for the nonzero fermion masses. This has been checked analytically and confirmed by 
the numerical calculation. We would like to stress at this point that this result is obtained 
with two independent widths T\y and Tz which violate the SU{2) gauge invariance. This 
finding seems to contradict the discussion of this issue in [0. The resulting violation of the 
high energy unitarity cancellations for e^er 4/ is suppressed by the factor Ty/Mw/s. A 
comparison of different gauge-boson width prescriptions performed in p[ and |12[ shows that 



our simple prescription is satisfactory as far as the experimental precision of LEP2 and future 
linacs is concerned. 



Using polarization vectors (pOf ) we can express the helicity amplitudes of reaction (|T9| ) corre- 
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spending to diagrams of Fig. |1| in terms of the functions defined in (|14|), (|T6j) and ( ]T8| ) 

M^ + M^ = F3 (m3, + g'-zJet, g^uB' + gz'Jsz^pt - pt2^P3 - Pu, "^3, 0, gw^w- ' 

Me + M7 = F3 (m3, 0, s-i^e^y- , - ^4.^72 - Pa^ "^4, 9^,de^ + 9-ide~ + 9"zd^~z^ ^4) , 

Mg + Mg = F3 (m5, + ^i^^^l, g^^,e~ + g^zl^'z^Pt - Pi2> - Pr2) "^5, 0, , , 

Mio = (M5,0,^iy£^y+,p^2 -Pe 5Pr2 -Pe >0>0'^2t'^z>^4) , (21) 

where we have used the shorthand notation Vi = fj(pj,Aj), i = 1,4,6, Uj = Uj{pj,\j), j = 
2, 3, 5 for the spinors which are defined according to @. The standard model couphngs of ( pOD 
and (^) are defined in terms of the electric charge e and the electroweak mixing parameter 
sin^ 9w- 

By comparing (^) with Fig. ^ one observes that the diagrams which differ only by the re- 
placement of the photon and Z propagators can be calculated simultaneously. This is one of 
the advantages of the presented method. 

Now the modulus squared of the spin averaged matrix element of the reaction ([T^) can be 
easily computed numerically. 

The matrix element of the bremsstrahlung reaction 

e+(pi) + e~{p2) u{p3) + d{pi) + /i"(p5) + ^^^(^6) + 7(^7), (22) 

where the particle four-momenta are indicated in parenthesis is calculated in the same way. 
The 71 Feynman diagrams of the process (^) can be obtained from those of Fig. [l| by attaching 
an external photon line to each charged particle as well as to the triple gauge boson vertex. 
We again neglect the Higgs boson contribution. 



In the soft photon limit, |p7| < u;, the matrix element of reaction ( |22D takes the simple 
factorized form 



P7|<w 



M^iPl,P2,P3,P4,P5,P6,P7, Ai, A2, A3, A4, A5, Ag, A7)| 
I Pi P2 , P3 Pi , P5 \ / ^ N /oQ^ 

- \9ii g-ii ^g-iu g-fd Vg^i \e^[pT,X-j) (23) 

V Pi- Pi P2- Pr P3 -Pr P4 ■ P7 Pb-PrJ 

X Mo(pi, P2, P3, P4, P5, P6, Ai, As, A3, A4, A5, Ae), 

where Mq is the matrix element of reaction (|1^) and the photon-fermion couplings are given 

by g^i = e, g^u = 2/3e and g^d = -e/3. 



4 The phase space integration 



The phase space integration is performed with the Monte Carlo integration routine VEGAS [|I5| . 
We integrate out the dependence on the azimuthal angle related to the rotational symmetry 
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with respect to the beam axis. This symmetry is satisfied as long as we do not consider 
transversely polarized initial beams. Thus, the number of integrations is reduced from 8 to 7 
and from 11 to 10 for reactions (0) and (P^, respectively. 

The 7 dimensional phase space element of the reaction (|19D is parametrized by 

7 A^/2(^s,S34,S56) ^^^^{s3A,mlml) X^/^{s56,ml,ml) 



d^Lips = {2tt 



8s 8S34 8S56 

X ds34ds56dcos6'd^73d^25, (24) 

where s = (pi +^2)^, -534 = {p^ +^4)^, s^q = (ps H-pe)^) ^ is an angle between the momenta 
Ps + P4 and the z axis of the cm. system which is directed along the positron momentum pi. 
df23 = dcos6'3d03 (dfis = dcosO^dcl)^ ) is the solid angle element of p3 (ps) in the respective 
cm. frame where Ps + P4 = (ps + pe = 0). 



The integration limits in the invariants S34 and S56 of (24) are given by 

(m3 + m4)^ < S34 < {^/s - - mef , (25) 
(m5 + maf < s^q < {Vs - v^)^ (26) 

and the spherical angles vary in the full range, i.e. 

-1 < cose<l, 

< fii < 47r, i = 3,5. (27) 



The 10 dimensional phase space element of the reaction (^2[) is parametrized in different ways 
dependent on whether we want to account for the peaking related to the photon emission from 
the initial or final state fermions. In order to deal with the initial state radiation peaking we 
parametrize the phase space by 

d^^Li s = f27rr^° — ^^^^^^^''^34,^56) A^/^(534, m^, m^) A^/^(g56, mg, m^) 
^ 2 8s' 8S34 8s56 

X d-E'7dfi7ds34ds56dcos6'34dfi3df25, (28) 

where s' = {pi + P2 — PiY- The photon variables, the energy Ej, and the solid angle Q7, are 
defined in the frame where pi + P2 — P? = 0. The polar angle ^34 of the momentum vector 
Ps + P4 with respect to the positron beam is defined in the same frame. The invariant masses 
S34, -556 and solid angles fl^, fl^ are defined as in 



The integration limits in the photon energy Ej and in the invariants S34, S56 of ( P5| ) read 

^cut <Er < (s - (m3 + 1714 + 7125 + me)^) /(2v^), (29) 
(m3 + m4)^<S34 < {y/s' - m5 - mef , (30) 
+ mef < S5Q < {^/s'-^/s^f, (31) 

where i^^cut is the minimum hard photon energy to be detected and s' = y/s{y/s — 2E-j). The 
spherical angles of (^8]) vary in the full range, as in 
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On the other hand, the radiation off the final state d quark is dealt with another parametriza- 
tion 



8s 8S56 

1 

X ■^dsa^jds^QdcosOsijdEsdEYdQ-sd^'sjdQ^, (32) 

o 

where the polar angle ^347 of the momentum P3 + P4 + P7 is defined in the c.m.s.; the energy 
£^3 and the spherical angle Q3 of the u quark, the photon energy Ej and the azimuthal angle 
037 between the u quark and the 7 are defined in the relative c.m.s. of u and d quarks and 
the photon; the spherical angle is defined in the c.m.s. of the /i~ and z/^. 

The integration limits are now specified by 

(m3 + 7714)^ < S347 < (v^-ms-me)^, (33) 
(m5 + m6)^<S56 < iVs-y/s^y, (34) 

r 1 1/2 

m3<^3 < [ml + \is,mlml{l + E,^jEr''))/{^s)\, (35) 



- 2(v/i-E3-|P3|) ' 



(36) 



where E'^^^ is the photon energy cut transformed to the c.m.s. of ud'^ and E'^'^^ = [m4+ 
A(s, m|, ml) I (4s)] is the maximum of the d quark energy. We have neglected E^cut on the left 
hand side of (|33D and E'^^^ in E^^^, which simplifies the corresponding analytic expressions. 
The correct phase space boundaries are then restored by checking the condition Ej > Ecnt in 
the c.m.s. numerically. The spherical angles of (|32|) vary again in the full range. 

The phase space parametrization convenient for the description of the photon radiation off 
the u quark or n~ is obtained from (|32|) by a permutation of the final state momenta. 

In order to improve the convergence of the phase space integration we perform the following 
mappings. The Breit-Wigner shape of the resonances is taken into account by the mapping 

sw = Twmw tan f ^^'^^ x + ipmin) + m^, (37) 



V iV, 



w 

where Nw is the normalization factor, Nw = ^wf^w/ ("^max — "^min), with ipnim = arctan(s{^™— 
m^r) / (Twfnw) and 4'rna.x = arctan(s,^'^'' — mly)/(Twmw)- The 1/t pole due to the the neu- 
trino exchange diagram (1) of Fig. ^ is mapped by transforming the polar angle of the virtual 
boson with respect to the positron beam 9w according to 

cosOw = (!-(! + Mrwn/f^w, (38) 
where Pw stands for for the velocity of the W'^ boson and = (1 + /3vi/)/(l — Pw)- 
The ^ l/Ej peaking of the bremsstrahlung photon spectrum is eliminated by the mapping 

E7 = Ef" (E^^''/Ef''y , (39) 
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where and E^^^ are the lower and upper hmit of the photon energy. The strong colhnear 
peaking behavior of the squared matrix element of reaction (0) corresponding to the radiation 
off the initial state positron ~ 1/(1 — /9 cos6'7) is eliminated by the mapping 



cos^7 = -^(l-(l + /3e)/r:) 



(40) 



with Te = (1 + Pe) /(I ~ Pe) sjid (3e = \J ~ 4:1711/ s being the velocity of the electron in the cm. 
system. Similarly, the collinear peaking related to the radiation off the initial state electron 
~ 1/(1 + /? cos 67) is dealt with the mapping 



cos^r = -^((l-/5e)r:-l). 

Pe 



(41) 



Finally the collinear and soft photon peaking corresponding to radiation off the final state 
fermion must be mapped away, e.g. for the d quark, l/(p4 ■ pr) ~ 1/(^*3 — -E3) is eliminated 
by the mapping 



Cs - {Cs - ms) 



C3 - 1713 



(42) 



where C3 = ^^347/2 + (m| — ml) / [2^/3347) . In ([37|-^2D , x denotes a random variable uniformly 
distributed in the interval [0, 1]. 



The phase space parametrizations (^) and (|32D together with the mappings (pTHi^) are 
implemented in a single multichannel Monte Carlo program in a way described in |T^. We 
use five different channels corresponding to the photon radiation off each charged particle with 
equal weights. 

In the soft photon limit, we can perform the integration over the photon phase space analyti- 
cally 



|d(j, 



1 



7l|p7|<a; 

P3 _ 
P3-P7 



(27r)3 

P4 



d^P7 ( Pi 

P7i<c^ 2i?7 y^^pi -pj 

2 



P2 



P2 -Pl 



g-yd 



Pi- Pi 



+ 9ji 



P5 



P5 -Pl, 



d(7o 



(43) 



The bremsstrahlung integrals J^^, which are defined by 



1 



d^Pr 



Pi ■ Pj 



27r)3 JiprKo; 2E7 {pi ■ pj){pj ■ p-j) 



(44) 



for i ^ i may be found in Sect. 7 of [0. For i = j we have 



, 2u 1 , 1 + A 
m — m 



(45) 



where A is the velocity of the radiating particle in the c.m.s. and rrij denotes a fictitious mass 
of the photon. 
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5 Numerical results 



We now present our results for the Born cross section e~^e ud^ and the corresponding 
hard bremsstrahlung process. We use the following physical parameters: the gauge boson 
masses and widths mw = 80.23 GeV, Tw = 2.085 GeV, mz = 91.1888 GeV, Tz = 2.4974 
GeV, the fermion masses: rUe = 0.51099906 MeV, = 105.658389 MeV, = 1777.05 MeV, 
TTiu = 5 MeV, 771^=10 MeV, ms=170 MeV, mc=1.3 GeV. The electroweak standard model 
couplings are parametrized by aw = 1/128.07 and by the electroweak mixing parameter 
sin^ 9w = 0.22591. The couplings of the bremsstrahlung photon are parametrized by a = 
1/137.0359895 which means in practice that we multiply the matrix element squared by the 
ratio a/aw- 

Our results have been thoroughly tested and checked against other calculations. The matrix 
elements have been checked against MADGRAPH I^Sf. Moreover, as we already mentioned in 
Sect. ^, we have checked the electromagnetic gauge invariance of the matrix element of the 
bremsstrahlung process both analytically and numerically. The phase space integrals have 
been checked against their asymptotic limits which have been obtained analytically. 



The total cross sections in the Born approximation are compared against EXCALIBUR ||16 
Tab. |1], with no cuts applied. 

Table 1: Born cross sections in pb (no cuts) 



m 



^cm GeV 


all 


< m 


162.5 


0.2685(0.4) 


0.2688(3) 


180.0 


0.6612(0.9) 


0.6616(7) 


189.0 


0.7037(1.0) 


0.7044(8) 


500.0 


0.2810(0.5) 


0.2817(5) 


1000.0 


0.1078(0.3) 


0.1079(2) 


2000.0 


0.03736(1.2) 


0.03748(8) 


10000.0 


0.002563(3) 


0.002578(16) 



In Tab. 0, we compare the Born cross sections ctq and the corresponding hard bremsstrahlung 
cross sections with the results of |]T2[ in the so called constant width scheme and with 
phase-space integration restricted by the "canonical" cuts. Let /, q, 7, and "beam" denote 
charged leptons, quarks, photons, and the beam (electrons or positrons), respectively, and 
0{i,j) the angles between the particles i and j in the cm. system. Furthermore, m{q,q') 
denotes the invariant mass of a quark pair qq'. The "canonical" cuts then read: 

^(/,beam)> 10°, 9(1,1') > 5°, 0(1, q) > 5°, 

^(7, beam) > 1°, ^(7,0 > 5°, 0{-f,q)> 5°, 

> 0.1 GeV, El > 1 GeV, E^ > 3 GeV, 

m{q, q') > 5 GeV . (46) 

Except for the additional angular cut between the charged leptons, which is irrelevant for the 
reactions considered here anyway, these cuts which exclude all collinear and infrared singular- 
ities coincide with those defined in fill. 
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Table 2: Cross sections^ in fb with the "canonical" cuts (^) 



^cm GeV 




^0 m 




a, H 


189.0 
500.0 
2000.0 
10000.0 


703.1(1) 
237.4(1) 
13.96(1) 
0.625(1) 


703.5(3) 
237.4(1) 
13.99(2) 
0.624(3) 


223.2(8) 
83.3(6) 
7.05(8) 
0.459(9) 


224.0(4) 
83.4(3) 
6.98(5) 
0.457(6) 



Another test, which is very sensitive to the treatment of the infrared and coUinear singulari- 
ties, is obtained by splitting the photon radiation cross section (^2]) into a soft and a hard part 



0"^ = (Ts + (Jh and checking the independence of the separation cut. The soft part includes the 



photons with Ey < uj and is given by 



The hard part then includes all photons with ener- 



gies Ej > LU. Furthermore, the total "inclusive" cross section is a = aQ + as + Ch. Since we have 
not yet included the infrared (IR) singular virtual one-loop corrections, as only exists when it 
is IR regularized in some way. Here we have chosen a small photon mass = 10~^ GeV. We 
demonstrate the cut (u) independence of the soft-hard splitting in Tab. ^ It is a measure of 
the numerical stability of our calculation as well as a test for the validity of the factorization 
into a radiation factor times the non-radiative cross section (|43|). 

Table 3: Cut-independence of = ag + (Th- The photon mass is = lO^^GeV^. 



^cm(GeV) 


uj (GeV) 


cTs (fb) 


ah (fb) 


as + ah (fb) 


189 


0.001 


202.6(2) 


1083(1) 


1285 




0.1 


712.2(5) 


572.8(6) 


1285 




1.0 


967.0(7) 


319.3(3) 


1286 


500 


0.001 


42.53(4) 


528.3(1.0) 


570.8 




0.1 


247.3(2) 


322.7(6) 


570.0 




1.0 


349.7(3) 


220.4(4) 


570.0 


2000 


0.1 


26.73(4) 


55.8(3) 


82.5 




1.0 


40.75(8) 


42.4(3) 


83.1 


10000 


0.1 


1.302(5) 


4.85(6) 


6.15 




1.0 


2.265(8) 


3.84(6) 


6.10 



Table 4: Mass dependence of cross sections (in fb; without cuts, except for E^ > 0.1 GeV) for 
different final states 



^cm GeV 


ao^udn'u^) 


aoicsn-u^) 


a^iudr Ut-) 


ay{udjj, u^) 


aj{csn'u^) 


a^iudr Ut) 


189.0 
360.0 
500.0 
2000.0 


704.1(4) 
422.0(2) 
281.0(2) 
37.33(4) 


703.8(4) 
421.8(2) 
280.9(2) 
37.32(4) 


703.5(4) 
421.5(2) 
281.0(2) 
37.32(4) 


573.4(4) 
448.5(4) 
322.8(4) 
56.48(27) 


525.2(4) 
418.4(4) 
302.0(4) 
53.19(25) 


522.6(4) 
414.1(4) 
298.1(3) 
52.67(13) 



^Here we adopt the parameters of [||: mw = 80.26 GeV, Tw = 2.05 GeV, mz = 91.1884 GeV, Tz = 2.46 
GeV, a = aw = 1/128.89 and sin^ 9w = 1 — rny^/m^. The fermion masses play no role in the presence of the 
cuts. 
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Finally, in Tab. ^ we illustrate the mass dependence of the cross sections of related channels. 
The replacements ud — ^ cs and — > r lead to comparable effects. While the Born cross 
sections remain practically unchanged the hard bremsstrahlung cross section, for the energy 
cut > 0.1 GeV, changes by about -9% (189 GeV) to about -6% (2 TeV). 



0.8 I — r 



0.7 - 



0.6 - 



0.5 - 



0.4 - 



0.3 - 



0.2 - 



0.1 - 



e+e udn 




{E^ > IGeV) 



J I I I I I I L 



_L 



200 400 600 800 1000 1200 1400 1600 1800 2000 



^ (GeV) 

Figure 2: The energy dependence of the total cross sections of reactions ( |19|) and (|2^) . 



The energy dependence of the total cross sections Uo and of reactions e~^e" udii~u^ and 
e^e~ udfi~u^'-f, respectively, is shown in Fig. |^. The hard bremsstrahlung cross section has 
been calculated with the photon energy cut of E^ = 1 GeV. 

At this point we would like to address the problem of SU (2) gauge-symmetry violation caused 
by introducing the constant widths F^y and Tz in a more quantitative way. For this purpose, in 
Fig. ^, we compare the e~^e~ —>■ udfi~u^ cross sections obtained when i) including all diagrams, 
ii) including only IVlV-diagrams and iii) assuming the creation of an on-shell IV^-pair and 
the subsequent decays W~ fJ'~i^fi and ud {e~^e~ — >■ W^W~ — >• udfi~Ufj,). Cases i) 

and ii) for the corresponding bremsstrahlung reaction with the photon energy cut E^ > 1 GeV 
are plotted in Fig. |^. Although the constant width prescription violates unitarity by spoiling 
the gauge cancellations in both cases i) and ii), the unitarity violation is much stronger in 
case ii) where we have neglected the non double-resonant diagrams. In case i) the effect is 
practically negligible, at least in the energy range presented in Fig. ^. This observation relies 



on the comparison with the results of [Q. Our results which were calculated in the linear 
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ft 



1 r-r 



0.1 



T" 



T 



T 



T 



e+e — > udfi 

a+ 



W+W- udfj.-u^ 




200 400 600 800 1000 1200 1400 1600 1800 2000 



(GeV) 

Figure 3: Lowest order cross sections obtained by including i) all diagrams, ii) only diagrams 
with 1^*1^* intermediate states and iii) production and decay of on-shell 1^-pairs (zero-width 
approximation) . 



b 




200 400 600 800 1000 1200 1400 1600 1800 2000 



(GeV) 

Figure 4: Bremsstrahlung cross sections for > 1 GeV obtained by including i) all diagrams, 
ii) only diagrams with W*W* intermediate states. 
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gauge agree within statistical errors with those of which were obtained in a nonhnear 
gauge in the so called complex-mass scheme that preserves the Ward identities. In Figs. |^ 
and ^ we observe that the double-resonant approximation provides a good approximation of 
the complete calculation for cm. energies from threshold up to about 250 GeV. However, it 
already deviates by 0.5% (more than 1% for reaction (p2|)) at ^/s = 500 GeV and by about 
2% at y/s = 1 TeV. On the other hand, we see that the zero width approximation, which is 
gauge invariant by definition, deviates from the complete tree level calculation by almost 18% 
at = 165 GeV, by 4.3% at = 200 GeV and by -4.5% at = 1 TeV. 



100 



10 



da 

I GeV J 



1 - 



0.1 - 



0.01 - 



0.001 




Vs = 500 GeV 



189 GeV 



50 



100 150 
(GeV) 



200 



250 



Figure 5: The photon spectra at ^/s = 189 GeV and ^/s = 500 GeV. 

The photon spectra at -^i = 189 GeV and = 500 GeV are shown in Fig. |. We see that 
the spectra are relatively soft, with a substantial fraction of events having photon energies of 
0{T\y). a bump of the 189 GeV spectrum at ~ 25 GeV reflects the ly-pair production 
threshold. We finally plot da/dml^^ at ^/s = 189 GeV as a function of the invariant mass of 
the ud pair 77134 in Fig. 
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T 




Figure 6: The differential cross section da/dm^^ at ^/s = 189 GeV versus the invariant mass 
of the ud pair 77134. The upper curve corresponds to the tree level reaction (^9[) and the lower 
curve to the bremsstrahlung reaction (^) with > 1 GeV. 

6 Conclusions and outlook 

We have presented an efficient method for calculating photon radiation cross sections for 
massive fermions. For the collinear region such finite mass calculations provide important 
tests for Monte Carlo generators which work with massless fermions. We have studied a 
complete signal plus background process with all possible real photon emission diagrams for 
the interesting channel e~^e~ ud^'v^. This channel is particularly suited for a detailed 
investigation of effects related to final state photon emission, since the muons appear well 
separated from photons in the detectors. In particular it seems to be interesting to study the 
influence of final state radiation on the W mass measurement via this channel. In addition at 
a high luminosity linear collider, like TESLA, one could study the quark mass effects due to 
the different quark flavor channels in e^e~ —>■ ^'v^^ + hadrons. Of particular interest would 
be a detailed investigation of the single top production channel e^e~ — > t6/i~z/^ which will be 
discussed in a forthcoming paper. 
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